Equation of State for Helium— 4 from Microphysics 



(N 
O 

o 

(N 

Or 
< 
00 



o 

a 

o3 



-a 
c 
o 

(N 
> 

in 



i 

C 

o 
o 



X 



Tim Gollisch 1,2 and Christof Wetterich 1 
1 Institut fur Theoretische Physik, Universitat Heidelberg, Philosophenweg 16, 69120 Heidelberg, Germany 
2 Innovationskolleg Theoretische Biologie, Humboldt- Universitat zu Berlin, Invalidenstr. 43, 10115 Berlin, Germany 

(April 8, 2002) 

We compute the free energy of helium-4 near the lambda transition based on an exact 
renormalization-group equation. An approximate solution permits the determination of univer- 
sal and nonuniversal thermodynamic properties starting from the microphysics of the two-particle 
interactions. The method does not suffer from infrared divergences. The critical chemical potential 
agrees with experiment. This supports a specific formulation of the functional integral that we have 
proposed recently. Our results for the equation of state reproduce the observed qualitative behavior. 
Despite certain quantitative shortcomings of our approximation, this demonstrates that ab initio 
calculations for collective phenomena become possible by modern renormalization-group methods. 
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I. INTRODUCTION 

When we regard a multiparticle system such as 
helium-4, we can choose between two possible points of 
view. We can zoom in on the details of the interactions, 
focusing on only a few particles, and in a first approach 
we might choose to study how the interactions of just two 
particles can be described. This is typically investigated 
by scattering experiments, and a quite precise picture of 
the two-particle interaction for helium-4 existsB On the 
other hand, we can also take a broader view and turn to 
the thermodynamic properties of many particles where 
properties .such as pressure, temperature, and superflu- 
idity arise. 

These two pictures are, of course, connected to each 
other as the thermodynamic properties are large scale 
effects of the (local) microscopic interactions. Yet for 
helium— 4 in the interesting region around the superfluid 
phase transition, a complete mathematical description, 
which makes this connection explicit, is still lacking, to 
our best knowledge. A general problem in the calcula- 
tion of thermodynamic properties is the occurrance of 
infrared divergences. Here, we present a novel approach, 
which is free of those divergences. It js based upon an 
exact renormalization-group equation.tj 

The goal of this work is threefold. First of all, we 
demonstrate how nonuniversal properties of thermody- 
namic systems can be calculated by means of an exact 
renormalization-group equation. Second, we present first 
results for the critical chemical potential and the equa- 
tion of state for helium-4 obtained merely from the mi- 
croscopic interactions. And third, we show the necessity 
of a shift of the chemical potential in the action, which 
arises from the mathematical manipulations of the func- 
tional integral and becomes important in the calculation 
of the phonon spectrum and the critical chemical poten- 
tial. We have reported on the appearance of this shift 
earlieiQ and now show its effect in application. 

The renormalization group is well suited for the treat- 
ment of the superfluid phase transition. While a pertur- 



bative approach is plagued by infrared divergences, the 
renormalization group naturally circumvents these diffi- 
culties. It also presents an intuitive connection between 
the microscopic and the macroscopic pictures. It starts 
with just the microscopic interactions. After integrating 
out the single-particle degrees of freedom consecutively 
from the effective ultraviolet (UV) cutoff at large mo- 
mentum down to zero momentum, a complete thermo- 
dynamic description of the system is obtained. 

This basic ideaQ is implemented in our approach by a 
functional differential equation that represents the lower- 
ing of an infrared cutoff. There is, in principle, no restric- 
tion on the relative interaction strength or the density of 
the system as in most perturbative approaches. Details 
of the method can be found elsewhere.! 

A similar approach to the calculation of nonuniver- 
sal properties for Bose condensates by renormalization- 
group methods was taken Jpy Bijlsma and Stoof, re- 
stricted to low density gases.cl Another recent application 
of renormalization group methods to superfluid helium 
combined with the exploitation of the underlying sym- 
metry and related Ward identity has successfully treated 
the infrared divergences—that occur in the computation 
of low-lying excitations £j 



II. THE ACTION FOR HELIUM-4 

We describe helium-4 by the many-body bosonic 
Hamiltonian operator 

H = (j^ _ a l a i + a li+q a l2- q v ( < l) a Q2a qi - 

1 ^ ' 31192,9 

(1) 

The mass m of the helium atom is 3.73 GeV, fi is the 
chemical potential, and v(g) the two-particle interaction 
potential. 
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Such a description is clearly not valid below the atomic 
scale, and we therefore assume a UV momentum cutoff 
proportional to the inverse atomic length scale given by 
the atomic diameter a: A = (27r)/er, q 2 < A 2 . (In all for- 
mulas, units are chosen such that h = c = ks = 1.) We 
choose to describe the interaction by a simple Lennard- 
Jones potential 



v(r) = 4e 



,12 



(2) 



with the maximum energy of attraction e = 10.22 K and 
a = 2.556 A. ThiSpfits the measured interaction poten- 
tial very accurately!!! For the definition of the interaction 
in the momentum domain, we identify the atomic di- 
ameter a with the scattering lengthEj and cut off the 
potential at small distances while fixing v(q = 0) = 
(4/ira) / (Vm) as required by low energy scattering theory. 
Here, V is the system's volume, which we will let even- 
tually approach infinity. The Fourier transform of the 
two-particle potential thus derived is shown in Fig. [l]. 
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v(q)+v(0)] 



(5) 



All momenta are restricted by the UV cutoff q 2 < A 2 . 

It is convenient to express the field as a sum over Mat- 
subara frequencies uj n — 2nn/ ' (3 
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FIG. 1. Fourier transform of the two-particle interaction 
for helium-4. 



This action differs in the appearance of Q from a 
functional-integral approach proposed earlier within a 
coherent-state ansatz, where a complete set of eigen- 
states of the annihilation operators is used for the copser- 
sion of the partition function to a functional integralO In 
contrast, Eq. (Q) is obtained by transforming the creation 
and annihilation operators in the Hamiltonian (Q) to lo- 
cation and momentum operators before the conversion 
to a functional integral. This allows us to choose an op- 
erator ordering that avoids uncontrollable contributions 
from next-to-leading-order terms in the discretization 
parameter. The correctness of our approach has beep, 
tested by numerical calculations in quantum mechanics.tJ 



The action <S[XjX*] f° r helium-4, which is a func- 
tional depending on the four-dimensional classical com- 
plex field Xj is defined via the grand canonical partition 
function 

Tr e- pH = J V X V X *e~ s[x ^ ] (3) 

with (3 = 1/T denoting the inverse temperature. As we 
have shown earlier,El great care has to be taken in con- 
structing the action, and for the Hamiltonian operator 
(|l|), we have obtained a perhaps unexpected shift £1 for 
the chemical potential /i in the action. One finds (up to 
a constant) 

s[x,x*] 

+ 2 '^2^qi+q^*q2-q( T ) v ( ( l)Xq2( T )Xqi w| (4) 

91:92,9 ' 



III. PHONON SPECTRUM 

From the action (fjj), we can already calculate the 
phonon velocity for a weakly coupled Bosc condensate 
at T = 0. Altough this can also be obtained by a Bo- 
goljubov transformation of the Hamiltonian operator (Q), 
we sketch the calculation as it shows how f2 cancels a cor- 
responding term evolving from the one-loop correction. 
In this way, fl acts as a counterterm though it has nat- 
urally emerged from the construction of the functional 
integral. 

The generating functional of the connected Green func- 
tions W[J, J*} is defined by 

e Wl= Jv X Vx*exp^-S[x,X*} 

F y ] {Jn,qXn,q + Jn,qXn,q) f ! ($) 
n,q ) 
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and the generating functional of the 1PI Green functions 
by the Legendre transform 



IV. FLOW EQUATION 



r[<P> <P*] = - hlW[J, J*} + ( J n,«^".« + J n,lfn,, 



(9) 



To calculate the phonon spectrum, we need the zeros 
of the inverse full propagator r^ 2 \ We restrict our- 
selves to the case of a static and homogeneous condensate 
¥n,q = tpdnflSq.o- Considering only low-momentum ex- 
citations, we can approximate the interaction potential 
v(q) by i>(0). The tree-level approximation yields 



r [<p,ip*] = s[<f,ip*] 

and therefore (for fixed n and q) 



(10) 
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I 5 2 S[ip,<p'} 



S<p—„ : — q6tp ni q 



iuj n + q v(0)tp* 2 
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v{0)ip 2 



-lU> r , 



where 



q = q 2 /{2m) -fi-Cl + 2v(0)\(p\ 2 . 
The one-loop correction is obtained as 

Fi[<p,<p* 



Itt inr' 2) 



o 



(11) 



(12) 



(13) 



with the trace including a sum of n and momenta. Ne- 
glecting terms of O (v(0) 2 ) , we obtain in the limit {3 — > oo 
the simple expression for the leading term, 



F (2) _ 
1 1 — 



(3tt 

o pn 



(14) 



This cancels the corresponding term in q appearing in 
Eq. (O) 



For + T\ z> , we can now compute the eigenvalues: 

„2 \ 



,(2) 
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2m 



/i + 2v(0) \tp\ J ± y/(iu; n ) 2 + v(0) 2 \v\ 



(15) 



The phonon energy spectrum is determined by the zeros, 
of A for complex oj n (with the first-order approximatior£3 
M = «(0)M 2 ): 

E(q) = iu} n (q) = \q\ 



v(0)\ip\ 2 +0(q 2 



(16) 



with the correct phone* velocity of a weakly coupled Bosc 
condensate at T = 0il3 



The physical situation of vanishing sources corresponds 
to an extremum of T since ST /Sip = J* . The value of T 
at the extremum T cq = r[<y9 eq , ip* q ] is directly related to 
the free energy F by 



F = TT 



eq 



fiN. 



(17) 



For a computation of thermodynamic quantities, we 
therefore aim at a computation of r cq as a function of 
T and /i. This involves a complicated functional integral 
over fluctuations (H|). We proceed by a stepwise solution 
by introducing a cutoff k for the fluctuations such that 
only fluctuations with momenta q 2 > k 2 are included 
in the functional integral. The variation of the effective 
action with k is described by an exact flow equation. 

For notational simplicity, we rescale the fields by ab- 
sorbing a factor of y/f3/ (2mV) into the field variable and 
define jl — 2m(fi + fi) and v(q) = 4m 2 Vv(q) / (3. For the 
action this yields 

S[x, X*} = V 2J Xn.q (2imuJ n + q 2 - p) Xn,q 
n, q 

,Y_ * * 

2 ' 4 ^n 1 ,q 1 +qXn 2 ,q2—q 

Til ,...,T14 

51 ,<22 -9 

(18) 

In the grand canonical partition function, we add 
source terms J and an infrared regulator Rk{n,q) and 
define the functional Wk[J,J*} as the generating func- 
tional of the connected Green functions in the presence 
of an infrared cutoff k, 



,W k [J,J*} 



J V X V X * exp j-S[ X ,X*] 

+ y ] \ Jn,qXn,q + Jn,qXn,q 



-VXn,qRk(n,q)x 



n,q , 



(19) 



For the infrared regulator, we choose 

Z k q 2 



Rk(n,q) 



exp(^)-l 



S n ,o- 



(20) 



Here, Z k is a wave-function renormalization, which will 
be defined below in Eq. (23). We note that this particu- 
lar regulator cuts off the momentum modes with q 2 < k 2 
only for the n = Matsubara frequency. The n ^ 
modes are not affected and contribute fully in the func- 
tional integral (f9h. For large enough T, the integration 



3 



of the n 7^ frequencies poses no problem since tem- 
perature acts as an infrared regulator for these modes. 
Nevertheless, a future extension of the infrared cutoff to 
the n^O Matsubara modes may be welcome for low T. 
The effective average action is defined as 



T k [tp, if*] = -W k + ^2 [ J n,«W 9 + J n,qV*n,q 



(21) 



As the infrared cutoff k is lowered, the evolution of Tk 
follows the exact renormalization group equatioro 



d k T k [<p*,ip] 



= iTr[(a fcJ R fc )(ri 2) +J R fe ) _1 ). (22! 



The functional differential equation ( |22| ) can only be 
solved approximately by truncating the most general 
functional form of the effective average action. For a 
truncation scheme, we choose a derivative expansion 



r fc 



^L 

P 



d 6 x 



U k {p) 



+ip*(r, x) ( 2m— - Z fe V 2 ) ip(r, x) 



--Y k (Vp)(Vp) 



(23) 



We consider an arbitrary dependence of the effective po- 
tential U k {p) on the invariant p(r,x) = tp*(T,x)(p(T,x). 
The derivative terms ( "kinetic terms" ) are multiplied by 
wave-function renormalizations Z k and Y k . The momen- 
tum and. field dependence of Z k and Y k are expected to 
be weako and therefore neglected for the computation of 
ground-state properties at this stage. For the calcula- 
tion of the excitation spectrum, these should turn out 
to be important and thus should be included in future 
calculations. 

We obtain evolution equations for the effective poten- 
tial and its derivatives by expanding the renormalization- 
group equation (|22| ) around a constant background field 
and taking the appropriate derivatives with respect to p. 
In particular, we will use the evolution equation for the 
first derivative of the effective potential U' k = dU k /dp, 
which is obtained as 



dtU k (p) 



(Z k q* + U' k +R k (q)y 



3U' k ' + 2U' k "p + Y k q 2 



(Z k q2 + Y k pq2 + U' k + 1U» 9 + R k (q)Y 



(24) 



Evolution equations for Z k and Y k are obtained by taking 
further appropriate expansions of Eq. (E2hJ§E-3 



The solution of U k for k — > is directly related to in- 
tensive thermodynamic quantities as pressure P, energy 
density e, or particle density n. Denoting by U oq (T,p) 
the value of U k =a at its minimum, one finds the relations 



P = ~TU e , 



dp, 



-T^+pn. (25) 



Therefore U cq (T,p) contains the information about the 
equation of state. Furthermore, the correlation length is 
encoded in the derivatives of U k= o(p) at the minimum. 
We therefore aim for a solution of the flow equation J24|) 
(or a similar equation for U k ) for k — > 0. 

This is obtained numerically by discretizing p and us- 
ing discrete approximations for the p derivatives of U k . 
As U k (p) remains fairly smooth throughout the evolu- 
tion, a relatively small number of sampling points suf- 
fice to give a good approximation of the shape of U k . 
The results presented in Sec. VI were obtained with 20 
sampling points in p space. In contrast to a Taylor ex- 
pansion of U k around its minimum, this retains more in- 
formation about the influence of regions of the effective 
potential away from the minimum. The resulting system 
of coupled ordinary differential equations was solved by 
an adaptive-stepsize Runge-Kutta integrator, where the 
momentum integrals were also evaluated numerically at 
each step. 

The generic evolution of U k (p) is discussed in detail 
elsewhere. □ In short, one observes that the minimum of 
U k (p) moves towards zero, but converges to a finite value 
Po for the superfluid phase with Uq(pq) = 0. Further- 
more, U k (p) is found to level out in the region p < po in 
accordance with the requirement that the true effective 
potential at k = must be a convex function of \ip\. For 
the nonsuperfluid phase, on the other hand, the minimum 
of Uq is at po = with Uq > 0. 

The finite number of sampling points is a possible 
source of numerical inaccuracies, especially in the ap- 
proach of the nonanalyticity at po in the superfluid phase 
and at the boundary p = in the nonsuperfluid phase. 
But since changing the number of sampling points does 
not substantially affect the results and since the criti- 
cal exponents are calculated with high accuracy as seen 
later, we do not expect these numerical sources to be 
major contributions to the errors of the results. 

The critical line separating the superfluid and non- 
superfluid phase corresponds to a "fixed-point" or 
"scaling" solution of the evolution equation in the 
dimensionless variables p = Z k k~ x p and u k {p) = 
k~ 3 U k (p). It is given by d k u k {p)\p = 0. Fig- 
ure H shows an example of the flow of the location 
k of the minimum of u k (p) for two different tem- 
peratures above and below the critical temperature. 
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FIG. 2. Flow of the location of the minimum of the dimen- 
sionless effective potential Ukift) close to the critical temper- 
ature (T c = 2.172 K, \T - T c \ w 1(T 7 K). One observes the 
approximate fixed point in the range ln(fc/A) « —4 to —10. 

We see that k stays close to the fixed-point solution 
for several orders of magnitude of k. Above the criti- 
cal temperature, it finally runs towards the symmetric 
phase where the minimum of the effective potential is at 
zero. Below the critical temperature, the minimum even- 
tually runs off in the other direction leading to a finite 
expectation value of the field and therefore to the case of 
spontaneously broken symmetry. 



V. DIMENSIONAL REDUCTION 

Since our cutoff acts only on the n = Matsubara 
frequency, the "initial value" of T k at k — A has to be 
computed as an integral over the frequencies. This 

procedure is known as dimensional reduction. It adds to 
the classical potential 



c/f {p) = -jip+\m? 



(26) 



a correction. The crudest approach would compute this 
in a one-loop approximation. In order to get a lit- 
tle more accurate, we choose, instead, a description in 
terms of a standard Schwinger-Dyson equation. This 
takes the leading mass renormalization into account. The 
Schwinger-Dyson equation is obtained by a one-loop cal- 
culation of the inverse full propagator of the Goldstone 
mode, with renormalized vertices given by Eq. (|23|). It 
yields a self-consistent equation, which we solve numer- 
ically at discretized points in p space for U'(p) while 
keeping the wave-function renormalizations and U"(p) 
fixed at their original values. This procedure is justi- 
fied a posteriori by the observation that the effect of the 
Schwinger-Dyson equation is nearly a constant shift of 
U'(p) for all p. 

The Schwinger-Dyson equation is obtained in the fol- 
lowing way. We decompose the field <p into its real and 
imaginary parts if = (tt + io~) / \/2 and write the trunca- 
tion as 



[tt, a] = J -j J d 3 xl U (p) + 2im ■ tt—o 

+ iz[(V 7 r) 2 + (Va) 2 ]+iy(Vp) 2 |. (27) 



The inverse full propagator of the Goldstone mode is then 
given by 



1 5T A [ir,a] 



= Z A q 2 + U' A (j>). 



(28) 



7T=0 



In order to calculate this propagator from the action Jig), 
we use the identity 



= / VnVa 



■ exp 



S [tt, a] 



E 

n . q 



—n,—q 



(29) 



where and J^ 7 ' are source terms for tt and a . This 
leads us to 







d d 



SS 



STTn.q [BjWdjM 



tO) \ P W[J (7 '\J < ' 7 '>] 
J -n,-q I e 



(30) 



We can now substitute the derivatives of W by the renor- 
malized vertices that we obtain from the truncation fl23| ) , 
e.g., by using the fact that is the inverse of the 
full propagator W^ 2 \ We denote the propagators of the 
Goldstone and the radial modes by P^ and P„: 

P, = Z k q 2 + U' k {p) + R k (n,q), 

P a = Z k q 2 + Y k pq 2 + U' k (p) +2pU"(p) +R k {n, q). (31) 

To account for the n dependence of the propagators, we 
distinguish the propagators of the zeroth Matsubara fre- 
quency (with R k (n, q) ^ 0) by writing P n and P a . 

To obtain the propagator of the Goldstone mode, we 
furthermore need to take the derivative of the resulting 

(tt) 

expression for Jl n _ q with respect to 7r_„._ g . We neglect 
two-loop graphs and explicitly evaluate the sums over 
Matsubara frequencies stemming from the 5S/6tt term 
in Eq. (^). This yields (for vanishing momentum and 
Matsubara frequency zero) 



1 



s 2 r A 



V 8-K n ^q8Tl^ n -q 
_ 1 8J- n -q 



7r=0,g=0,ra=0 



VStt- 



TT=0,q=0,n=0 
A 



-£ + 5(0)p+JL J dqq 2 (P^ 



5 



1 



v(0) + 2i(q)] P a + v(0)P« ~ 4: P U'l(p)v(p) 

1 



dq q 



21 L^oth^VKK 



yj Ptt Per PttPq 

[«(0) + 2%)] P a + v(0)P* - ApU'l{p)v{p) 



(32) 



By equating the right-hand sides of Eqs. (28) and (32) 



we obtain the Schwinger-Dyson equation, which we solve 
numerically for U' A (p) with U'j[(p) — v(0) on a discrete 
set of p values. This yields the initial condition for the 
evolution of U' k (p) . 

A similar procedure could be used for a determination 
of Z\(p) and Ya( P ) from the Schwinger-Dyson equation 
at nonzero q 2 . We omit here the fluctuation effects for 
the wave-function renormalization and use the "classi- 
cal" values Za = 1, Y\ = 0. We believe that this sim- 
plification gives a sizeable contribution to the error in 
the computation of nonuniversal quantities presented in 



Sec. VI 



We are aware that the use of the Schwinger-Dyson 
equation becomes problematic for low T. In this region, 
one would prefer to use an infrared regulator that drops 
the Kronecker delta in Eq. (ppf). As a result, also the 
contributions of the n ^ Matsubara frequencies to the 
Schwinger-Dyson integral would be suppressed by the 
cutoff A, and all modes with q 2 < A 2 could then be 
dealt with by the renormalization-group procedure (in 
contrast to the present version where this holds only for 
n = 0). One would thus expect a more reliable treat- 
ment. Modifying the cutoff ( pfj| ) by dropping 8 n $ results 
again in a simple form of the evolution equation since the 
Matsubara sum can be evaluated explicitly. We find 



d k U k (p) 



d 3 q 



(2tt) 3 
(3 P. + Pc 



[d k R k {q)} 



(3 



&m^/P~p; 



:C0th-f-v/p~~pT. 

4m 



(33) 



A numerical solution of this evolution equation, however, 
shows instabilities, which we have not managed to over- 
come so far. 



findings of Maynard,0 where a value of around — 7.4 K 
has been extrapolated from the measurements of the ve- 
locity of the fourth sound in the superfluid phase. We 
emphasize that the Q shift in Eq. (Q) is essential, as oth- 
erwise we would have obtained a value of about 12 K 
higher. This supports again our ansatz for the action for 
helium-4. 

Critical exponents. By exploring the region around 
the critical temperature, we can compute the critical ex- 
ponents v, (3, and 7 of the helium-4 system. They de- 
scribe the divergence of the correlation length above and 
below the critical temperature (v), the growing of the 
order parameter, i.e., the expectation value of the field, 
below the critical temperature (/?), and the divergence 
of the susceptibility above and below the critical tem- 
perature (7). We present our results as a demonstration 
that the flow equations can indeed be followed directly 
from microphysics to macrophysics without encountering 
infrared problems. No rcsummations of series or other 
technical tricks as used by other analytical methods are 
needed. 

It is well established that helium-4 at the lambda tran- 
sition is described by the 0{2) universality class. The 
critical exponents can therefore be compared to previ- 
ous high-precision estimates for this universality class. 
The excellent correspondence with other, more demand- 
ing calculations and experimental values shown in Ta- 
ble I demonstrates that our truncation is reliable for 
the universal quantities. The scaling requirement that 
v and 7, respectively, obtain the same values above 
and below the critical temperature is also well fulfilled. 

TABLE I. Critical exponents for the superfluid phase tran- 
sition of helium-4 and comparison with values for the 0(2) 
universality class. The slight differences in the values obtained 
for v and 7 above and below the critical temperature reflect 
numerical inaccuracies. 
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0.675 a 
0.667 b 



0.359 b 



0.6695 c 
0.671 d 
0.672 6 

0.3455 c 
0.3485 d 



0.67095(13) 
0.6705(6) 
0.6708(4) 



VI. RESULTS 

From the solution of the evolution equation, we can 
calculate the quantities of interest for the phase transi- 
tion to superfluidity. 

Critical chemical potential. For a constant temper- 
ature of 2.172 K, which corresponds to the experimental 
critical temperature under vapor-pressure conditions, we 
can tune the chemical potential to obtain the fixed-point 
solution. We thereby find a critical chemical potential 
p c = —6.7 K in good agreement with the experimental 



1.305 3 
1.298 b 



1.316 c 
1.315 d 
1.33 e 



a For T > To. 
b For T < T c . 

c From summed perturbation series at six-loop order. 
d From a fifth-order e expansion. 
c From lattice calculations. 

Pressure. Since ST /Sep* = J, the minimum of the 
effective potential Uk=o{po) corresponds to the physical 
case of vanishing external sources for the stationary he- 



G 



lium system. The pressure P is given by the value of the 
effective potential at the minimum po, P = —TUk=o{po)- 
The minimum value U k (po(k)) can be followed by an evo- 
lution equation, which is obtained from Eq. (J22j). It is 
given by 



d k (u k (p (k)) 
d 3 q 



(27T) 



dkRk{q) 



Z k q 2 
l 



Rk(q) 



(34) 



Z k q? + Y k p (k)q2 + 2U^po(k))Po(k) + Rk(q\ 
in the spontaneously broken regime (po(k) 7^ 0) and by 
d 3 q d k R k (q) 



d k (U k (0)) 



(2irf Z k q2 + U> k (0) + R k (q) 



(35) 



in the symmetric regime (po(k) = 0). The starting value 
is the classical minimum J7a(po(A)) = — p 2 /[2v(0)]. 
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FIG. 3. Temperature dependence of the pressure for three 
different values of the chemical potential fi calculated from 
the evolution equation (|24|). (The value of /1 = —6.74 K cor- 
responds to the chemical potential at the critical temperature 
under vapor pressure as obtained from the fixed-point solu- 
tion.) 

We find that the solution of the flow equations does not 
lead to a reliable estimate of the absolute value of the 
pressure. The main reason for this is a strong depen- 
dence on the UV cutoff, which is not well controlled by 
our method. For example, the initial value of C/a(po(A)) 
should also receive a correction from higher Matsubara 
frequencies, which we do not calculate in our Schwinger- 
Dyson approach. This correction will depend notably on 
the microphysical details of the interaction and the choice 
of the UV cutoff. Fortunately, this is not the case for the 
temperature dependence of the pressure, which is mainly 
governed by the momentum range q 2 ~ 2mT. We thus 
resort to normalizing the pressure at the critical temper- 
ature, where p is given by the fixed-point condition, and 
only calculate the T dependence of P. Figure || shows the 
dependence of the pressure on the temperature for three 



different values of the chemical potential. Since we need 
to specify the chemical potential in our formalism when 
we want to calculate the order parameter and the density 
at vapor pressure, we determine p{P, T) from this result 
and use the experimentally known vapor-pressure curve 
for P(T). 

Order parameter. While the critical exponent (3 de- 
scribes the growing of the order parameter in the 0(2) 
universality class, we can also calculate the order param- 
eter for the helium system explicitly. It is given by the 
expectation value of the field (p. Therefore, po(k = 0) 
gives the magnitude squared of the expectation value of 
the field. In Fig. g we show the results for the tem- 
perature dependence of po = pa(k = 0). The super- 
fluid density is proportional to pq. The proportionality 
constant, however, depends on the wave-function renor- 
malization, which is not accurately determined so far. 
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2.1 2.2 
Temperature [K] 
FIG. 4. Magnitude squared of the order parameter of the 
superfluid phase obtained from the evolution equation (pi|). 
The chemical potential is adapted so that the calculated pres- 
sure yields the vapor pressure corresponding to the tempera- 
ture. 

Density. Finally, we compute the particle density as 
1 /• 



V 2V 



yTr(r^ )~\ (36) 



The inverse full propagator is block diagonal in n and q 
with 



r ( k \n,q) = 



s 2 r k j v , v *] \ 



(37) 



and contains for k — the information about the phonon 
energy spectrum for given p and T and arbitrary strength 
of the interaction. The trace in Eq. ([36]) includes a sum- 
mation over momenta and Matsubara frequencies. We 
calculate the eigenvalues of from the truncation (|2 
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and explicitly sum over the Matsubara frequencies. This 
yields (in the rescaled field variable) 

(38) 

where P n and P a again denote the inverse propagators 
of the Goldstone and the radial mode as in Eq. (|3l|), but 
with Rk=o(n, q) = 0. The masses Uq and Uq + 2U^'p are 
taken at po, so that in the spontaneously broken phase 
we have Uq(pq) — and in the symmetric phase po = 
and Pyr — P a . 

As in the calculation of the pressure, the calculations 
of the density do not give reliable results for the ab- 
solute value of the density, which strongly depends on 
the UV cutoff and is expected to get a sizeable con- 
tribution from three-particle interactions. We observe 
that the momentum integral in Eq. (|38| ) is dominated by 
high momenta q 2 « A 2 and is therefore strongly affected 
by the choice of the cutoff. Furthermore, this implies 
that the precise momentum dependence of the propaga- 
tors P T and P a could also have a substantial influence 
on the value of the density. The temperature depen- 
dence of the density, on the other hand, is dominated by 
q 2 2mT. We thus allow for an overall shift of the den- 
sity and only calculate its temperature dependence. The 
results are shown in Fig. pj. They reproduce the qualita- 
tive picture of the density close to the phase transition. 
The peak at T c is clearly visible as well as the slightly 
negative thermal-expansion coefficient for T < T c . 



P a should be replaced by momentum-dependent func- 
tions Z(q 2 ,T), Y(q 2 ,T). We are aware that the simpli- 
fications Z(q 2 ,T) = Z\ — 1 and Y(q 2 ,T) = Y A = are 
less appropriate for the temperature dependence of the 
density, for which the integral is dominated by q 2 ~ 2mT. 
In the vicinity of the phase transition, the contribution to 
dn/dT from the fluctuation integral ( pq ) scales roughly 
oc Z~ 3 / 2 . Replacing Z\ = 1 by Zk T , k^= 2mT leads to 
a less pronounced temperature dependence if Z^ T > 1. 
Furthermore, for T near T c , critical fluctuations influ- 
ence the behavior of Z and Y at low q 2 . We believe that 
our crude treatment of the wave-function renormaliza- 
tion constitutes the most important source of error in 
determining the temperature dependence of the density. 
We nevertheless see that we obtain qualitatively correct 
results, which are stable against shifting the UV cutoff. 

Since the temperature dependence of the density is 
not a universal critical quantity, an understanding of the 
momentum and temperature dependence of Z(q 2 ,T) is 
mandatory for a quantitative prediction. This is illus- 
trated by Fig. [| where we show the temperature depen- 
dence of a momentum-independent Z(T) that would be 
needed for a reproduction of the experimental data. (The 
latter may be viewed as an effective T-dependent mean 
value of Z(q 2 , T), which yields the same q 2 integral ( pSSj ) 
as the true Z(q 2 ,T).) This curve seems reasonable, in- 
cluding the modest feature near T c . Using the dotted 
curve of Fig. |^ as an approximation to these Z(T) val- 
ues still leads to considerable deviations from the ob- 
served density in the immediate vicinity of T c , demon- 
strating the high precision to which Z has to be known. 
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FIG. 5. Temperature dependence of the density. The den- 
sity is calculated for different temperatures around the criti- 
cal point (filled boxes) and compared to measurements (solid 
line). The chemical potential is adapted so that the calcu- 
lated pressure yields the vapor pressure corresponding to the 
temperature. 
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FIG. 6. Temperature dependence of the wave-function 
renormalization Z (filled circles) that would be needed in 
Eq. ( p8| ) to yield the observed density. The dotted line dis- 
plays a simple exponential fit to the data. 



The fact that our calculations deviate substantially 
above the critical temperature is due to the crude treat- 
ment of the wave-function renormalization Z. In the 
present approach, we have used Z(q 2 ,T) = Z\ = 1 and 
Y(q 2 ,T) — Y\ — 0. For a more reliable treatment, the 
wave-function renormalizations Zk=o and lfe=o in P v and 



Another question why more precision on the wave- 
function renormalizations will be needed concerns the 
calculation of the excitation spectrum and the sound ve- 
locity. This requires an understanding of the q 2 depen- 
dence of the inverse full propagator T^ 2 ) and thus strongly 
depends on Z(q 2 ,T) and Y(q 2 ,T). 
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VII. CONCLUSION 

We have computed thermodynamic properties of a 
strongly interacting system near a second-order phase 
transition from "first principles" . Starting from the mi- 
crophysical interactions, we have determined the macro- 
scopic thermodynamic potential for helium-4 near the 
superfiuid phase transition. Our method is based on 
an exact renormalization group equation. The critical 
exponents obtained by this approach agree well with 
more elaborate state-of-the-art calculations and mea- 
surements. Our most robust nonuniversal result, the crit- 
ical chemical potential, is also in good agreement with 
experiment. It supports our computation of the micro- 
physical action by which helium-4 is described in the 
language of functional integrals. This action contains an 
additive shift of the chemical potential that has not been 
considered so far. 

For quantitatively convincing results of the equation of 
state, one needs to improve the treatment of the UV cut- 
off as well as additional effects of higher Matsubara fre- 
quencies beyond those that we have taken into account in 
the Schwinger-Dyson formalism. A more accurate treat- 
ment of the wave-function renormalization is required as 
well. Qualitatively, though, we reproduce important fea- 
tures of the equation of state such as the density peak at 
the phase transition. 

We emphasize that we treat here a temperature range 
where collective effects and critical behavior are crucial 
and where many previous methods are plagued by severe 
infrared problems. Despite certain quantitative short- 
comings, this work may be viewed as a demonstration 
that ab initio computations of the equation of state be- 
come possible with modern renormalization-group tech- 
niques. 
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